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Dark energy with the equation of state w(z) rapidly evolving from the dustlike (w ~ at z ~ 1) 
to the phantomlike (—1-2 < w < — 1 at z ~ 0) has been recently proposed as the best fit for 
the supernovae la data. Assuming that a dark energy component with an arbitrary scalar-field 
Lagrangian p(ip, V pip) dominates in the flat Priedmann universe, we analyze the possibility of a 
dynamical transition from the states (ip,ip) with w > — 1 to those with w < — 1 or vice versa. We 
have found that generally such transitions are physically implausible because they are either realized 
by a discrete set of trajectories in the phase space or are unstable with respect to the cosmological 
perturbations. This conclusion is confirmed by a comparison of the analytic results with numerical 
solutions obtained for simple models. Without the assumption of the dark energy domination, 
this result still holds for a certain class of dark energy Lagrangians, in particular, for Lagrangians 
quadratic in V^ip. The result is insensitive to topology of the Priedmann universe as well. 



I. INTRODUCTION 

One of the greatest challenges in modern cosmology 
is understanding the nature of the observed late-time 
acceleration of the universe. The present acceleration 
expansion seems to be an experimental fact, now that 
data from supernovae type la 0, 0], corroborated later 
by those from the cosmic microwave background @], have 
been recently confirmed by the observations of the largest 
relaxed galaxy clusters Q. Although the observations are 
in good agreement with the simplest explanation given by 
a cosmological constant A of order (1CT 3 eV) 4 , the myste- 
rious origin of this tiny number, which is about 120 orders 
smaller than the naive expectations, gives rise to the idea 
of a dynamical nature of this energy. Possible dynami- 
cal explanations of this phenomenon are given in various 
frameworks. One of them is known as quintessence (see 
e.g. (2^1 and other references from the review j24[). In 
this framework the equation of state p = we is such that 
W > — 1. Another proposal is the phantom scalar fields 
(see e.g. which possess the super-negative equation 
of state w < —1, due to the "wrong" sign before the 
kinetic term in the Lagrangian. Alternatively, there is 
a more general possibility under the name k— essence 
0- EH HH which is an effective scalar-field theory de- 
scribed by a Lagrangian with a nonlinear kinetic term. 
For this model, the equation of state w is not constrained 
to be larger or smaller than —1. Allowing the dark energy 
to be dynamical provides an opportunity to study the so- 
called coincidence problem which asks why dark energy 
domination begins just at the epoch when sentient beings 
are able to observe it. The main advantage of k— essence 
is its abilit y to solve this problem in a generic way (for 
details see PI|]), whereas the first two models require a 
fine-tuning of parameters. 

Without imposing the prior constraint w > — 1, 
the observations seem to favor the dark energy with 
the present equation of state parameter w < — 1 (see 
e.g. Ref. El El El US])- Moreover, recently it was argued 
(see Ref. P-Lal an d other constraints on w(z) obtained in 



Refs. |1 Oil III III III 113 ) that the dark energy with the 
equation of state parameter w{z) rapidly evolving from 
the dustlike w ~ at high redshift z ~ 1, to phantomlike 
— 1.2 < w < — 1 at present z ~ 0, provides the best fit for 
the supernovae la data and their combinations with other 
currently available data from the measurements of cos- 
mic microwave background radiation (CMBR) and from 
2dF Galaxy Redshift Survey (2dFGRS). 

Matter with w < — 1 violates the dominant energy con- 
dition which is a sufficient condition of the conservation 
theorem j2fj|. Therefore for such models one cannot guar- 
antee the stability of vacuum on the classical level. The 
instability can reveal itself at the quantum level as well. 
In fact, it was shown that the phantom scalar fields are 
quantum-mechanically unstable with respect to decay of 
the vacuum into gravitons and phantom particles with 
negative energy [ll], LL2]- Assuming that the phantom 
dark energy is an effective theory allows one to escape this 
problem through the appropriate fine-tuning of a cutoff 
parameter. If the dark energy could dynamically change 
its equation of state from a phantomlike one to that with 
W > — 1, then this transition would prevent the undesir- 
able particle production without such a fine-tuning. Here 
it is worth mentioning that quantum effects on a locally 
de Sitter background could lead to the effective parame- 
ter w < -1 (see Ref. |3flLl4oj). 

Another fundamental physical issue where this tran- 
sition could play an important role is the cosmological 
singularity problem. If w < —1 in an expanding Fried- 
mann universe, then the positive energy density of such 
phantom matter generally becomes infinite in finite time, 
overcoming all other forms of matter and, hence, leads to 
the late-time singularity called the "big rip" [14] . The 
transition under consideration could naturally prevent 
this late-time singularity. Here it is worthwhile to men- 
tion that for certain potentials and initial conditions the 
phantom scalar fields can escape this singularity by evolv- 
ing to a late-time asymptotic which is the de Sitter so- 
lution with w = — 1 \?i4 \. Moreover, it was argued 
that the quantum effects can prevent the developing of 
the "big-rip" singularity as well [3S| . 
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On the other hand, to avoid the big crunch singularity, 
which arises in various pre-big bang and cyclic scenarios 
(see e.g. |2Ml2^l23), one assumes that the universe can 
bounce instead of collapsing to the singularity. The ex- 
istence of a nonsingular bouncing solution in a flat (or 
open) Friedmann universe (k ^ +1) requires the viola- 
tion of the null energy condition (e + p > 0) during the 
bounce j^jj. If the energy density s is constrained to be 
positive, then it follows that w < — 1 is the necessary 
condition for the bounce. But the energy density of such 
phantom matter would rapidly decrease during the col- 
lapse and therefore only the transition from w > — 1 to 
w < — 1 just before the bounce could explain the non- 
singular bouncing without a fine-tuning in initial energy 
densities of phantom and other forms of matter present 
in the universe. 

It is worth noting as well that for regimes where the 
equation of state of the k— essence field is greater than 
— 1 it is possible to find a quintessence model which gives 
the same cosmological evolution but behaves differently 
with respect to cosmological perturbations Hence, 
it is interesting whether this equivalence can be broken 
dynamically. 

In this paper we consider the cosmological dynamics of 
a A;— essence field p, described by a general Lagrangian p 
which is a local function of <p and V M y>. The Lagrangian 
depends only on ip and a scalar quantity, 



x = iv^vv 



(LI) 



First of all, we determine the properties of a general La- 
grangian p(<p,X), which are necessary for the smooth 
transition of the dark energy from the equation of state 
w(ip,X) > —1 to w{ip, X) < —1 or vice versa. The tran- 
sition obviously happens if the system passes through the 
boundaries of the domains in the space (ip, X), defined 
by these inequalities. In most of the paper, we assume 
that the dark energy dominates in a spatially flat Fried- 
mann universe. The main question is whether trajecto- 
ries connecting these domains on the phase space (ip, tp) 
exist and are stable with respect to cosmological pertur- 
bations. In the case of the phase curves which do not 
violate the stability conditions, we study their asymp- 
totic behavior in the neighborhood of the points where 
the transition could occur. To proceed with this analysis, 
we linearize the equation of motion in the neighborhood 
of these points and then use the results of the qualita- 
tive theory of differential equations. For the dark energy 
models described by Lagrangians linear in X, we perform 
this investigation beyond the linear approximation. For 
this class of Lagrangians, we illustrate the outcome of 
our analysis by numerically obtained phase curves. Fi- 
nally, we generalize the results to the cases of spatially 
not-flat Friedmann universes filled with a mixture or the 
dark energy and other forms of matter. 



II, GENERAL FRAMEWORK 

Assuming the dominance of the dark energy, we neglect 
all other forms of matter and consider a single scalar field 
<p interacting with gravity. After all, we will see that 
the results can be easily extended to the models with 
additional forms of matter. The action of the model reads 



in our units (M p = H 



1, where M p is the reduced 



Planck mass M p = (8ttG)-^ 2 = 1.72 x 10 18 GeV ) as 
follows: 



S — Sn ~\~ S<£ 



d 4 x, 



(II.1) 



where R is the Ricci scalar and p(<p, X) is the Lagrangian 
density for the scalar field. This kind of action may 
describe a fundamental scalar field or be a low-energy 
effective action. In principle, the Lagrangian density 
p((p, X) can be non-linear on X. For example, in string 
and supergravity theories nonlinear kinetic terms appear 
generically in the effective action describing moduli and 
massless degrees of freedom due to higher order gravita- 
tional corrections to the Einstein-Hilbert action (l7tllfll]. 
The "matter" energy-momentum tensor reads 



T, 



6S^ 



5g 



{11.2) 



= p.x(p, X)Vfj,pV v (p - p(<p, X)g llu . 



Here a comma denotes a partial derivative with respect 
to X . The last equation shows that, if V v p is timelike 
(i.e. X > 0), the energy- momentum tensor is equivalent 
to that of a perfect fluid, 

Tfu, = {e + p)UfJJ v - pg^, (II.3) 

with pressure p(<p, X), energy density 

e{ip, X) = 2Xp iX (cp, X) - p(<p, X), (IL4) 

and four velocity 



(11.5) 



The equation of motion for the scalar field can be ob- 
tained either as a consequence of the energy-momentum 
tensor conservation V M T^ = or directly from the ex- 
tremal principle SS v /5p = 0: 

p, x U g p + p.xx (V^V^) WV "tp + e, v = 0, (II.6) 

where D g = .g^V^V" and V M denotes the covariant 
derivative. For this fluid we can define the equation of 
state parameter w as usual: 



P 

w — —. 

e 



(II.7) 



There is increasing evidence that the total energy den- 
sity of the universe is equal to the critical value, and 
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hence in the most part of the paper we will consider a 
flat Friedmann universe. In the end, we shall show that 
the results are also applicable in the cases of closed and 
open universes. Thus, the background line element reads 

ds 2 = g^dx^dx" = dt 2 - a 2 (t) dx 2 . (II.8) 

The Einstein equations can be written for our back- 
ground in the familiar form: 

£ = ~(e + 3p), (II.9) 

H 2 = |, (11.10) 

where H = a /a is the Hubble parameter and a dot 
denotes derivative with respect to the physical time t. 
These equations also imply a continuity equation: 

E= -3H(e + p). (11.11) 

In general, whenever n / any two of these three last 
equations imply the third one (by compatible initial con- 
ditions). Usually it is easier to work with the second 
and the third equations (these are the Friedmann equa- 
tions). Note that, from Eq. ijll.lOJl . e was constrained to 
be non-negative. 

Because of the homogeneity and isotropy of the back- 
ground, we get X = \ip 2 and p^ = <pp : x so the en- 
ergy density looks as the energy in usual ID classical 
mechanics 

e(tp,ip) = <pp,<j>-P- (11.12) 

Expressing H from the first Friedmann equation i|II.10Jl . 
we can rewrite Eq. (Ill . Gil in the case of the homogeneous 
and isotropic flat background as follows: 

<pe <x + <pp,xV3s + e, v = 0. (11.13) 

So far as a(t) ^ 0, all of the information about the dy- 
namics of gravity and scalar field is contained in the equa- 
tion written above. In accordance with our initial simpli- 
fication the dark energy should dominate in the universe; 
therefore we assume throughout the paper that e > 0. 

Following we introduce the effective sound speed 
of the perturbations, 




(11.14) 



Then the equation of motion takes the form 

ip + ipc 2 s V3e + — = 0- (11.15) 

In most of this paper, we shall assume that the solutions 
ip(t) and Lagrangians p(tp, X) have enough continuous 
derivatives. So, for example, (p(t) will be mostly consid- 
ered as being at least of the class C 2 : p>(t), ip(t), <p(t) 
are continuous. 



III. POSSIBLE MECHANISMS OF THE 
TRANSITION 

There are two possibilities for the evolution of dark 
energy from w > — 1 to a phantom dark energy with 
w < — 1 (or vise versa). These are a continuous transi- 
tion, in which the dark energy evolves through points 
where w = — 1, and a discontinuous transition occurring 
through points where e = 0, provided that the pressure 
p is finite. Since by assumption the dark energy is the 
dominating source of gravitation, we cannot have e = 
and therefore it is sufficient to consider only continuous 
transitions. 

Further, throughout the paper we will usually suppose 
that for the dynamical models under consideration there 
exist solutions tp(t) and corresponding to them moments 
of time t c such that 

w[<p(t c ),X(t e )]=-l. (IH.l) 

Henceforth the index c denotes a physical quantity taken 
at t c ; i.e., p c = <p(t c ), e c = e{t c ), etc. 

The parameter w can be expressed with the help of 
Eq. iTT.4ll in the following form, more convenient for a 
study of continuous transitions: 

2X 

w = -l+ px- (III.2) 

£ 

Since X = ^ip 2 > and e > 0, we find that w < — 1 
corresponds to p t x < 0, whereas w > — 1 implies 
p t x > 0. In accordance with our notation, the equation 
of state parameter w takes the value —1 at the points 
= (cp c , ip c ), where either X = or p t x = 0. Because 
of these equations, the points ^ c generally form curves 
7(A) and isolated points in the phase space (<p, tp) of the 
dynamical system given by Eq. AI.13L The curves 7(A) 
may intersect. 

For our purposes, it is convenient to divide the set of 
points 'J'c into three disjunct subsets: 

A) The tp— axis of the phase plot (tp, tp), i.e., ip c = 0. 

B) The points where p,x{^c) = but ip c 7^ and 

£,x(* c )^0. 

C) The points where p,x(^c) = and s,x(^c) = but 

Pc^O. 

Further in this section we will study the dynamics of the 
scalar field p in the neighborhoods of ^ c separately for 
these cases. If the system evolves from the states (<p, <p) 
where w > — 1 to the states with w < — 1 (or vise versa), 
the function p t x changes sign. 

It is worth noting that, if the scalar dark energy were 
equivalent to an "isentropic" fluid for which the pressure 
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p is a function only of e, then the possibility of evolving 
through the points e c where w(e c ) — — 1 could be eas- 
ily ruled out. Indeed, in that case we could rewrite the 
continuity equation ijll.lli only in terms of e: 



(HI 3) 



3e [e+p(e)], 



so that the system of Einstein equations l|II.9j) and ijll.lfljl 
could be reduced to Eq. (IIII.3I) and the values of energy 
density e c would be fixed points of this equation. A dy- 
namical transition through a fixed point is clearly impos- 
sible. 

An example of the dark energy which seems to be 
equivalent to the "isentropic" fluid is the simple model 
described by the Lagrangian p = p(X) depending only 
on X. Let us further assume that there are some values 
X c , where w(X c ) = — 1, If Eq. HI.4II is solvable with 
respect to X in the neighborhoods of these points X c , 
then one can find X(e) and therefore the pressure is a 
function only of energy density p(e) = p[X(e)]. Thus 
the system is equivalent to the "isentropic" fluid, X c are 
fixed points, and the transition through w = — 1 is im- 
possible. It remains to consider the conditions on the 
function e(X) under which Eq. (Ill . 41) is solvable with re- 
spect to X. From the theorem about the inverse function, 
Eq. l|II.4Jl is solvable with respect to X if 

e, x {X e ) = [2Xp tXX (X c )+ P}X (X c )} + 0. (IIL4) 

One can see directly from the equation of motion l|II.13Jl 
and condition llIII.4ll that X c are fixed point solutions. In 
fact, as it was shown in Ref. (2^1, there generally exists 
the solution X{t) = X c and moreover it is an attractor in 
an expanding Friedmann universe. Thus, the transition 
is generally forbidden for systems described by purely 
kinetic Lagrangians p(X). 

In the general case when p = p(tp,X), the pressure 
cannot be expressed only in terms of e, since tp and X 
are independent. 



A. Transition at points X c = 

Here we will analyze the possibility of the transition 
in the case (A) cp c = Q. Namely, we are going to study 
the properties of the solutions ip(t) in the neighborhood 
of the line tp c = 0. Differentiating the equation of state 
parameter with respect to the time, we have 



2X 



2X 



-P,x 



-P,x 



2X 



(III.5) 



At the points under consideration we have w c = be- 
cause X c = and, respectively, X c — tp c tp c — 0. More- 
over, the time derivatives in the second and third sum- 
mands vanish at these points as well due to the continuity 
equation ijll.llll and the formula 



Let us differentiate the w once more with respect to the 
time. The only term which survives from the formula 
ljIII.5l) at the points X c — is the first term. Hence, we 
have 



2X 



-P,x 



2(p 2 



-P,x 



(111.7) 



Using the equation of motion (III. 13|i . we can express (p 
through the p and its derivatives 



(j'C 



(III.. 



p'x = tp (p, v x + PP.xx) ■ 



(HI 6) 



As follows from Eq. l|II.4Jl . e tX (t c ) = P,x(t c ) at the time 
t c when the system crosses X c = 0. Provided £, v (t c ) ^ 
and p t x(t c ) ^ 0, we infer from Eq. 11111.71) and ijlll 8|) that 
the equation of state parameter w(t) has either a mini- 
mum or a maximum at the point t c . Thus, the transition 
is impossible in this case. 

If £,<p{t c ) — and £, x (t c ) ^ 0, then it follows from 
relation (IIII.8|) that tp(t c ) = 0. Therefore the considered 
solution ip(t) for which w [(p(t c )(p(t c )] = — 1 is a fixed 
point solution (p(t) = cp c = const and the transition is 
impossible. Since e > 0, we see that this fixed point is 
obviously the de Sitter solution. 

If not only X c = but also e.x(t c ) = P,x(t c ) = 0, 
then s. v (t c ) — and it follows from the formula (111.-41) 
that p^(t c ) — 0. Moreover, the equation of motion 
llll.l.'tll is not solved with respect to the highest deriva- 
tives (namely, with respect to (p ) and therefore does not 
necessarily have a unique solution. It happens because 
the point (<p c , 0) on the phase plot (tp, (p) does not deter- 
mine the ip via the equation of motion ljll.131 . It is clear 
that, in this case, the pointlike (on the phase plot) solu- 
tion tp(t) = ip c = const is a solution, but not necessary a 
unique one. 

Below, we will give a more general consideration of the 
geometry of phase curves in the neighborhood of the <p> 
axis. The phase flows are directed from right to left for 
the lower part of the phase plot ip < and from left to 
right for the upper part (p > 0, see Fig^ Therefore the 
system can pass the <p> axis only if the point of intersection 
is a turning point (curves 1 and 6 on Fig . Otherwise 
the crossing is a fixed point (or a singularity) . If there is 
a smooth phase curve on which the system pass through 
the (p axis, then in a sufficiently small neighborhood of 
the turning point we have (tp — <p c ) ~ l p> 2n ~ X n where 
n > 1. Restricted to this curve, the function Xp t x(f, X) 
depends only on X and in the absence of a branching 
point the sign of this function above and below the <£ 
axis is the same. Then it follows from the formula (IIII.2j) 
that the system cannot change the sign of (w + 1) while 
crossing the ip axis. 

If a smooth phase curve does not cross but touches 
the ip axis at a point ip c (see Fig. trajectories 3 and 
4), then the following asymptotic holds: (p ~ (tp — tp c ) 2n , 
where n > 1. Let us find the time needed for the system 



5 




C. Transition at points \& c : p,x(^c) = 0, 
e,x(*c) = 0, X c 

As we have already mentioned at the beginning of this 
section, the points ^ c generally form the curves in the 
phase space (tp,tp). The subclass of the points which 
we are going to consider in this subsection, is generally a 
collection of the isolated points given by the solutions of 
the system, 



Figure 1: Possible phase curves in the neighborhood of the 
tp— axis. Only on the curves 1 and 6, the system crosses the 
tp axis. Curves 2, 3, 4 and 7 have an attractor as a shared 
point with the tp axis, whereas curves 5 and 8 have a repulsor. 
These attractors and repulsors can be fixed-point solutions or 
singularities. 



to reach the tangent point (<p c ,0) in this case. We have 

rVc dip r- dp 



p,x(p, <p) = 0, e tX (p, <p) = 0. 



(HI 11) 



t : 



(<P-<fc) 



2n ' 



(III.9) 



where ipi n is a starting point on the phase curve. The 
last integral is obviously divergent. Therefore the system 
cannot reach the tangent point in a finite time 

Finally, we come to the conclusion that in the frame- 
work under consideration it is impossible to build a model 
with the desirable transition through the points X c = 0. 



B. Transition at points * c : p,x(*c) =0, 
e,x(*c) 7^0, X c 

In the neighborhood of a point \& c , at which the con- 
dition E t x(^c) 7^ holds, one can find a function tp c (tp c ): 
P,x [tfc, tPciVc)] = 0. This follows from the theorem 
about the implicit function. One would anticipate that 
on the phase curves intersecting the curve tp c (tp c ) the 
state of the dark energy changes to the phantom one 
(or vise versa). Let us express p : x from Eq. IIIII.2I) and 
substitute it into formula (III. 141) for the sound speed of 
perturbations: 



(w + l)e 
2Xe 



x 



(111.10) 



For the stability with respect to the general metric and 
matter perturbations the condition c 2 s > is necessary 
(see Indeed the increment of instability is inversely 

proportional to the wavelength of the perturbations, and 
hence the background models for that C 2 . < are vio- 
lently unstable and do not have any physical significance. 
Because of the continuity of e,x, there exists a neighbor- 
hood of the point ^ c where e t x ^ 0. Therefore, from the 
above expression for the sound speed IIIII.IOII it follows 
that if (w + 1) change a sign then c 2 s should change a sign 
as well. If this is the case, then the trajectories, realiz- 
ing the transition, violate the stability condition (? s > 0. 
Therefore the model of the transition is not realistic. 



Only for specific models, the solutions of this system are 
not isolated points. An example when these solutions 
form a line is considered in section Hvl Usually the phase 
curves passing through the isolated points build a set of 
the zero measure. Therefore it is physically implausi- 
ble to observe the processes realized on these solutions. 
The only reason to study the behaviour of the system 
around these points is their singular character. The point 
is that the equation of motion HI, 1311 is not solved with 
respect to the highest derivatives at this points. In such 
points there can be more than one phase curve passing 
through each point. Moreover, the set of solutions <p(i), 
which pass through ^ c with different tp, could have a 
non-zero measure. On the other hand, the equation of 
motion does not necessarily have a solution tp(t) such 
that ((p(t c ),(p(t c )) = at some moment of time t c , or 
there exists the desirable solution ip{t) but it does not 
possess the second derivative with respect to time at the 
point ^c. Below we will analyze the behavior of the phase 
curves in the neighborhoods of the points ^ c . 

can be rewritten as a 



The equation of motion ill.U 
system of two differential equations of the first order: 



d(p 
~dt 

dtp 

~dt 



-ip- 



, P,X 



3e 



£,x' 



(111.12) 



tp. 



The phase curves of this dynamical system are given by 
the following differential equation: 



dtp 
dp 



tpp^xVSs + 
tpe,x 



(111.13) 



This equation follows from the system JIII.12II and there- 
fore all phase curves corresponding to the integral curves 
of system IIIII.12I) are integral curves of the differential 
equation l|III.13Jl . But the reverse statement is false, so 
each integral curve of Eq. I|III.13J| does not necessarily 
correspond to a solution tp(t) of the equation of motion 
(or of the system l|III. 121) 1. In the neighborhoods of the 
points where e t x ^ 0, it is convenient to introduce a new 
auxiliary time variable r defined by 



dt = e.xdr. 



(III. 14) 
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The system ilIII,12Jl is equivalent to the r system: 



dip 
dtp 



(111.15) 



The auxiliary time variable r change the direction if e t x 
change the sign. Note that the system IIIII.15I) always 
possesses the same phase curves as the equation of motion 
HTTT.12L 

In the case under consideration we have ip c ^ and 
from the formula 



dtp (p 
dip ip 



(111.16) 



we infer that dtp/d(p(t c ) should be finite, if <p(t c ) is finite. 
As one can see from the equation determining the phase 
curves i|III,13Jl . in order to obtain a finite dip/ dip (t c ) it 
is necessary that at least e iV (^ c ) = 0. In the case if 
£,v{^c) 7^ 0, the solution ip(i) does not possess the sec- 
ond t derivative at the point t c . Usually this can be seen 
as unphysical situation. But nevertheless this does not 
necessarily lead to the unphysical incontinuity in the ob- 
served quantities e, p, H, and ip, (p, One may probably 
face problems with the stability of such solutions, but 
let us first of all investigate the behavior of the phase 



curves in the case e )V (W c ) 7^ 0. From Eq. ijIII,13}) . we 
obtain dip /dip = at $ c . Further we can parameterize 
the phase curve as ip — ip(ip) and bring the equation for 
phase curves l|III.13Jl to the form 



dip 

-7T = e, x {y,<p)F{ip,ip), 
dip 



where we denote 



F(ip,ip) = -- 



■•P 



<PP,XV<l£ + £ 



(111.18) 



If, as we have assumed, e )¥ ,(^ c ) ^ and e^c) 7^ 0; then 
F(tp, ip) is differentiable in the neighborhood of the point 
and F(^ c ) — —ip/e >ip . For the second ip derivative at 
the point H! c , one obtains 



d^p 
d 2 ip 



£,xx- 



(111.19) 



That is why the point ^ c is a minimum or a maximum for 
the function ip(ip). In this case f c is such an exceptional 
point on the phase plot, where the solution ip(t) cannot 
have continuous ip(t) and the phase curve terminates (see 
points £ in Figs. Eland EJ. This happens because the di- 
rection of the phase flow is preserved in the neighborhood 
of * c . If e y xx(^c) = 0, then one can find the third ip 
derivative of <p(<p) at the point 



d\ 
d 3 ip 



sxxx- 



(111.20) 



In this case there can exist a continuous solution ip(t) 
such that (ip(t c ), tp(t c )) = <3> c at some moment of time t c 
and the only bad thing happening in this point is that 
ip(t c ) does not exist. Let us now investigate what hap- 
pens with the equation of state at this point of time. Dif- 
ferentiating both sides of the definition l|II.7Jl of w yields 



<P t 



(111.21) 



where we have used the equation of motion IIII.13j) at 
the point t c and the definition l|II.14JI of (? s . Applying 
the l'Hopital rule for the c 2 s (t c ) 



lim p x Is x , we find 

t-*t e 



that cj(t c ) = 0. Moreover, using the l'Hopital rule 



at the point ~$> c , one can find 
Thus, if p iV 7^ the transition 



for the derivative of c s 
that del /dip ~ P,v/ £ ,v- 
could occur but it changes a sign of the sound speed c 2 . 
Therefore, if the stability criteria are applicable to this 
case, then the transition leads to instability. 

The necessary condition for the existence of ip during 
the transition is 

e,„(*c) = 2X c p,xA^c) - pA^c) = 0. (111.22) 

This condition drastically reduces the set of the points 
\& c , where the transition is possible. Namely, they are 
the critical points of the function e(ip, X) and, on the 
other hand, they are the fixed points of the auxiliary 
t system l|III.15Jl . These fixed points are additional 
to the fixed points of the system 



(111.17) = 0, e iV (cp,0) 



mT2| defined by 
0, and e t x{ip,0) ^ 0. From now 
on, we will consider only those points where the 
condition l|III.22Jl holds. From the relation l!III.22j) . it 
follows that if p yip (ip c ,X c ) — then p,x v {<p c , X c ) = 0. 
Otherwise X c — 0, and as we have already seen the 
transition cannot happen via the points X c — 0. Note 
that if p,(p(^t) — 0, then the points are common 
critical points of the pressure p(ip,(p), energy density 
e(ip,ip), and p^xi^P,^)- From condition 1III.22tl follows 
that points \E'+ are singular points of Eq. 13|) . In 
such points there can be more than one phase curve 
passing through this point. Moreover, as we have already 
mentioned, the set of solutions tp(t), which pass through 
\P+ with different ip, could have a nonzero measure. 
For example, if were a nodal point (see Fig. HJ, 
there would be a continuous amount of the solutions 
passing through this point and therefore there would be 
a continuous amount of solutions on which the transition 
could occur. 

Let us investigate the type of the singular points 
This will tell us about the amount of the solutions ip{t) on 
which the transition is possible and their stability. For 
this analysis, one can use the technique described, for 
example, in 2l|, and consider the integral curves of the 
equation l|IIL13JI . Here we proceed with this analysis in 
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a more convenient way, namely, using the auxiliary t sys- 
tem . It is convenient because for this system the 
singular points are usual fixed points. As we have 
already mentioned, both systems have the same phase 
curves and therefore the analysis to perform is also ap- 
plicable to the phase curves of the system l|III.12Jl . The 
only thing we should not forget is the difference in the di- 
rections of the phase flows of these systems. Ife(5'+) 7^ 0, 
then one can linearize the right-hand side of the r system 
in the neighborhood of a point V E'+: (ip^ + S(p, y>+ + Sip). 
The linearized r system l|III.15Jl is 

-^-V = AV, (111.23) 

UT 

where we denote 

and elements of the matrix A are given by the formulas 

a = £,js><P, (111.25) 

b = 2Xe, xx , 

c = - (3H(pp iXv > + £ w ) , 

.9 = -£,x v ¥>, 

where all quantities are calculated at Here we have 
used the Friedmann equation ijll.lflll . If is an iso- 
lated fixed point of the r system flll.151 (or equivalently 
the singular point of system l|III.12Jl ). then the following 
condition holds 

detA = ag-6c^0. (111,26) 

The type of the fixed point depends on the eigenvalues 
A of the matrix A (for details see, for example, (13]). In 
the case under consideration a — —g and therefore we 
have 

A 2 = be + a 2 = - det A. (111.27) 

If be + a > 0, then eigenvalues A are real and of the 
opposite signs. In accordance with the classification of 
the singular points, \E'+ is a saddle point (see Fig. [5J. 
Therefore the transition is absolutely unstable; there are 
only two solutions cp(t) on which the transition is allowed 
to occur. 

If be + a 2 < then A are pure imaginary. Here the 
situation is a little bit more complicated: In accordance 
with [2]]], this fixed point of a nonlinear system can be 
either a focus or a center. In these cases, as one can see 
from Fig. there no solutions ip(t) passing through the 
point Therefore, from now on we will consider only 
the first case - real A of the opposite signs. 

It is convenient to rewrite the expression for A into 
a simpler form. Differentiating the continuity equation 
(IrTTTl yields 

e c = -3H c p c . (111.28) 




Figure 2: Phase curves in the neighborhood of the singular 
point are plotted for the case of the real A. At the points 
£, the solutions <f>(t) do not exist. These points together with 
form the curve F on which e,x(r) = 0. 




Figure 3: Phase curves in the neighborhood of the singular 
point are plotted for the case of the pure imaginary A. 
Here we assume that the singular point is a focus. 



Remember that the index c denotes quantities taken at t c 
or in this subsection at Differentiating the pressure 
p as a composite function, we have p — p^p + p,xX. 
Assuming that (p c is finite, we obtain that p c = P^Pc- 
Thus the formula for e c is 

E C = SHcp^pc (111.29) 

Using the condition JIII.22II and the last equation, we 
bring the element c of the matrix A to the following 
form: 

c=-< w + ^- (HI.30) 

This relation allows one to rewrite the formula l|III.27Jl 
as follows: 

A 2 = 2X C (( £ ^) 2 _ e c ^ xx ) + e c e< xx . (111.31) 
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Figure 4: If A had the eigenvalues Ai = A2, then the singu- 
lar point HFJ would be a nodal point and there would be a 
continuous set of trajectories passing through it. To illustrate 
this we plot here the phase curves in the particular case of a 
degenerate nodal point. The form of the equation of motion 
HI. 13(1 excludes such types of singular points and therefore 
prevents the possibility of such transitions. 

Here it is interesting to note that the expression 

e %ip £C xx ~ ( e ?x v ) from the previous formula is the de- 
terminant of the quadratic form arising in the Taylor set 
of e in the neighborhood of the critical point ^+ . If this 
determinant is positive, then the function s(ip, X) has 
either a minimum or a maximum at point Other- 
wise there is either one curve of constant energy density 
e(curve) = with a singular turning point at 

in the neighborhood of ore two intersecting at ^+ 
curves of constant e (see |2£|) ■ On the other hand, dif- 
ferentiating £ as a composite function we find 

e c = 2X C (s c >vv + 2e c >Xv Cp c + e c >xx ip^) . (111.32) 

Substituting this relation into the previous formula 
for A yields 

A 2 = 2X C (e° xx <p e + £ c Xv f . (111.33) 

This formula provides the relation between ip c at the mo- 
ment of transition and A. Note that A depends on (p c only 
in the case when e xx (t c ) = 2X c p xxx (t c ) 7^ 0. More- 
over, comparing formulas l|III.29Jl and l|III.32Jl for e' c , one 
can obtain the equation on (p c : 

^xxVcvl + 2E e tXv ip c ipe + + ZH cP c v = 0. (111,34) 

This equation is solvable in real numbers if the discrim- 
inant is positive. As one can prove, the discriminant is 
exactly the 4A 2 and therefore positive if we consider the 
saddle point. The same can be seen from relation l|III.33Jl 
as well. 

Let us denote the positive and negative eigenvalues and 
the corresponding eigenvectors of A as A + , A_ = — A + 
and a + , a_, respectively. 



If b ^ ( or equivalently s c xx 7^ 0), then the eigen- 
vectors can be chosen as a + = (1,(A + — a)/b) and 
a = (1, — (A + + a)/b). Therefore the separatrices form- 
ing the saddle are 

A+ — a \+ + a 

dip, = dp and dip_ = dp. 

b b 

The general solution for the phase curves in the neigh- 
borhood of V E'+ is 

fc- ^+ ~ a ? \ f c ■ ^++ a c \ 

[dip — dip J I dip A — dtp \ = const. 

(111.35) 

If b = and additionally a > 0, then we have A + = 
a = £,X(p<j>) an d one can choose the eigenvectors as 
a + = (l,c/2a) and a = (0,1) . For a negative a, one 
can obtain the eigenvectors and eigenvalues by changing 
A+ <-* A_ and a + <-> a_. The separatrices then read 

dp = and Sip — —5p. 

2a 

Thus, similarly to the previous case, the phase curves are 
given by 

Sip [Sip — —S(p) = const. (III. 36) 

V 2a I 

As we have already mentioned from the formula 
dip/dip= tp/v, follows that, at the points where the 
phase curves are parallel to the <p axis and where ^0, 
the second t derivative of the field does not exist. If 
we look at the equations l|III.35|l and l|III.36Jl providing 
the phase curves in the neighborhood of 5'+, then we 
find that in the first case (when e c xx 7^ 0) the phase 
curves lying on the right- and left- hand sides of both 
separatrices should have a point £ where they are par- 
allel to the ip axis (see Fig. |2j . Therefore each of these 
phase curves consists of two solutions of the equation 
of motion l|II.13Jl and the exceptional point £ where the 
solution ip(t) does not exist. The same statement holds 
in the case of the pure imaginary A (see Fig. UjJ. This 
behavior is not forbidden because, as one can easily 
prove, the exceptional points £ lie exactly on the curve 
r on which £,x(r) = and the equation of motion is 
not solved with respect to the highest derivatives. Note 
that we have already assumed A / 0, and from this 
condition it follows that s c xx and £ c X(p cannot vanish 
simultaneously. Therefore in the neighborhood of the 
point there exists an implicit function (p(tp) (or 
<p(ip)) and its plot gives the above-mentioned curve T on 
which £,x(r) = 0. In the case £° xx = 0, the separatrix 
Sep = locally coincides with T, and therefore, this 
integral curve of Eq. l!III.13jl does not correspond to any 
solution ip(t) of the equation of motion. Nevertheless, 
in virtue of the existence theorem, all phase curves 
obtained in the neighborhood of the separatrix Sip = Q 
correspond to the solutions of the equation of motion. 
Moreover, if e c xx = and P c Xip 7^ then the curves 
on which p x = and e x = locally coincide with 
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each other and with the curve Sip = 0. The only phase 
curve intersecting the curve Sip = at is the second 
separatrix Sip = cSip/2a. Thus the only solution ip(t) 
on which the transition happens in the neighborhood of 
corresponds to the separatrix Sip = cSip/2a. This 
can also be seen from the equation l|III.34Jl which has 
only one root ip c in this case. In the Sec. ijlVJl we will 
illustrate this with a numerical example (see Fig.^J. 

It is worthwhile to discuss cases which fall out from 
the consideration made above. We have assumed that 
A 7^ and therefore *f>+ is an isolated singular point 
of Eq. The most natural possibilities to drop 

out this condition are s c Xlp = and either e c xx = or 



ip c 3H c p 



<P<P 



= 0. In the first case ^fj is a critical 



point not only of the function e but of the function 
e t x as well. This can be obtained either for a very 
special kind of function p namely, such that p,x = 0, 
p,xx = 0, p,xxx = 0, and [p - 4X 2 p. X x] = at 
or imposing the condition that the point is a critical 
point not only of p but also of the functions p,x and 
P,xx- p,x = 0, p, v = 0, p,xx = 0, p,x v = 0, p,xxj> = 0, 
and finally p,xxx — at In the second case 

is a common critical point for the functions e and £ ><p , 
In terms of p, this condition is as follows: p,x = 0, 
P,<P = 0, p^Xip = 0, p t xx = 0, p iVV = 0, p,xx v = 0, 
and finally p.xipip = at . Thus, the point V E'+ is a 
common critical point of p, p t x and p t xip- Of course, the 
analysis performed above does not work in the case if 
the function p(ip, X) does not have a sufficient amount of 
derivatives. It is clear that all these cases are not general. 

Let us sum up the results obtained in this section. In 
the general case of linearizable functions e t x, £,<p, and 
p t X, the considered transitions either occur through the 
points \E , +, where p.x = 0, s t x = 0, e. v = 0, and 

<P [{ £2 x v ~ £, vv £ : xx) V - 3He^xxP, v ] > 

or lead to an unacceptable instability with respect to 
the cosmological perturbations of the background. The 
points ^+ are critical points of the energy density and 
are the singular points of the equation of motion of the 
field ip as well. These singular points are saddle points 
and the transition is realized by the repulsive separatrix 
solutions, which form the saddle. Therefore the measure 
of these solutions is zero in the set of trajectories and the 
dynamical transitions from the states where w > — 1 to 
w < — 1 or vise versa are physically implausible. 



IV. LAGRANGIANS LINEAR IN X 

The simplest class of models, for that one could antic- 
ipate the existence of dynamical transitions, is the dark 



energy described by Lagrangians p(<p, X) linear in X: 

p{ip, X)=KX-\v=\ (K(p)V^ip - V(ip)) . 

(IV.l) 

In the isotropic and homogeneous Friedmann universe, 
the Lagrangian is then 



V(ip)) 



(IV.2) 



For these models, we always have <? s = 1 and therefore, as 
follows from our analysis, the transitions could occur only 
through the points where £ : x — 0- The energy density 
for this model is 



e(ip,ip) = ±(K(ip)ip 2 + V(ip)) 



(IV.3) 



If one takes K(ip) = 1, then the Lagrangian l|IV.2|l is 
the usual Lagrangian density for scalar field with a self- 
interaction. If we take K(ip) = —1, then we obtain the 
so-called "Phantom field" from and [[J. The case 
K(ip) > corresponds to w > —1, whereas K(ip) < 
corresponds to ic < —1. The equation of motion ill. 13tl 
takes in our case the following form: 



(IV.4) 



ipK + pKVSe + e v = 0. 



While the equation determining the phase curves ljIII.1 
reads in this particular case 



^ + + -^—e tV = 0. 
dp ipK 



(IV.5) 



If K(ip) is a sign-preserving function, one can redefine 
field ip: y/\K(ip)\dip = d(j> (see also Ref. ^f|). The equa- 
tion of motion for the new field <j> can be obtained from 
Eq. ljrV4Jl . through the formal substitutions ip — ► <f>, 
V(ip) -> V((f>) ee V(ip(<p)), and K(<p) -> ±1, where the 
upper sign corresponds to a positive K{ip) and the lower 
one to a negative K{ip). After these substitutions, the 
equation of motion ilIV.4Jl looks more conventionally 



( 



±& + V((f>) 



, l(dV(4>)\ 
2\ d<t> J 



0. (IV.6) 



Moreover, this equations is easier to dial with, because 
one can visualize the dynamic determined by it, as ID 
classical mechanics of a point particle in a potential 
±V(4>)/2 with a little bit unusual friction force. If we 
were able to solve the equation of motion ijlVfiJl for all 
possible V(<p) and initial data, we could solve the problem 
of cosmological evolution for all linear in X Lagrangians 
with sign-preserving K{ip). 

If the function K(<p) is not sign-preserving, then at 
first sight it seems that the dark energy, described by 
such a Lagrangian, can realize the desirable transition. 
The function K(ip) generally can change the sign in two 
ways: In the continuous one, then the function K(ip) 
takes the value zero for some values of field ip or in a 
discontinuous polelike way. 
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A. Linearizable K{tp) 

Without loss of generality, one can assume that 
K(Q) = K c = 0, K(ip) < for the negative values of 
cp and K(p) > for p > 0. The line p = on the 
phase plot (p, <p) we will call the "critical" line for the 
given class of Lagrangians. The phantom states ((p, ip) 
of the scalar field lie on the left-hand side, while the 
usual states with w > — 1 are on the right-hand side of 
the critical line. If there exists a solution tp(t) whose 
phase curve passes through the "critical" line, then the 
dark energy can change the sign of (w + 1) during the 
cosmological evolution. From now on, we will investigate 
the behavior of the phase curves of the system in the 
neighborhood of the critical line. 

First of all, it is worth considering the functions K(cp) 
such that K' c > (here we have denoted K jtp (0) = K' c ), 
because in this case we can directly apply the outcome of 
our previous analysis made in the Sec. IIII Al Condition 
(IIII. 221 is for the linear in X Lagrangians as follows: 



p 2 c K' c + V' c = 0. 



(IV.7) 



As we have already assumed K' c > 0, therefore, if V' c > 0, 
then, as follows from condition IIIV.7II . there are no twice 
differentiable solutions <p(t) whose phase curves would 
intersect or touch the critical line. Further (see formula 
l|IV.15ll and below) we will show that, for the linear in X 
Lagrangians, condition llIH.221 (or in our case condition 
IIIV.7II is necessary not only for the existence of the second 
t derivative ip at the point of intersection with the critical 
line but for the existence of a solution ip{t) at this point 
as well. Thus, we come to the conclusion that, if Vj. > 0, 
then two regions <p < 0, and <p > on the phase plot are 
not connected by any phase curves and accordingly the 
dark energy does not change the sign of (w + 1) during 
the cosmological evolution. 

In the case V' c < 0, we can solve Eq. Ill V. 71 with respect 
to tp c : 



(IV. 



The phase curves, lying in the neighborhoods of the sin- 
gular points Hf£ = (0, u±), are to obtain from the relation 
ijTTTm which gives: 



where 



p \ tp — u± -^-tp j = con*! . 



nrT V'K'J-Vl'K' 
A+ = -3H C + ■ c ■ ' 



2u ± (K' c ) 



(IV.9) 



(IV.10) 



For each singular point (0, u±), there is a corresponding 
solution p±(t) whose phase curve is the separatrix 



A± 

p ± =u± + —ip, 



(IV.ll) 



which intersects the critical line. These phase curves 
correspond to the const — in the right-hand side 
of Eq. IIIV.9|) . Another curve, which corresponds to 
const = is ip — 0. As we have already mentioned at the 
end of the previous subsection, this curve does not cor- 
respond to any solutions p(t) of the equation of motion 

(frOf . 

Considering the phase flow in the neighborhoods of 
\f + = (0, u±) (see Fig. EJ, we infer that the separatrices 
p± are repulsors immediately before they intersect the 
critical line and attractors after the crossing. Hence, the 
measure of the initial conditions (ip, ip) leading to the 
transition to phantom field (or vice versa) is zero. In this 
sense the dark energy cannot change the sign of K(<p) (or 
equivalently the sign of (w + 1)) during the cosmological 
evolution. 

The typical behavior of the phase curves in the neigh- 
borhood of the singular points (0, u±), for the models 
under consideration (K' c > 0, V' c < 0), is shown in Fig. El 
Here, as an example, we have plotted the phase curves 
obtained numerically for a toy model with the Lagrangian 

density p = \<pp 2 — \ {{ip — l) 2 + |J . For this model we 
have u± = ±1, A + = — §, and A_ = — 5. 

Let us now consider such potentials V(p) that V£ = 0. 
If K{p) is a differentiable function, then, in the case un- 
der consideration, the equation of motion i|IV.4t obvi- 
ously has a fixed-point solution ip(t) = but this solu- 
tion is not necessarily the unique one. When V' c = and 

the 



K' c > 0, then, as follows from the condition l|III.22JI 
only value <p, where a phase curve could have coinciding 
points with the critical line, is ip = 0. From the analy- 
sis made in Sec. IIII Al we have already learned that the 
transition is impossible in this case. Nevertheless it is 
worth to showing explicitly how the phase curves look 
at this case. Taking into consideration only the leading 
order in the numerator and denominator of Eq. l|rV.5j) 
and assuming that V" ^ 0, we obtain 



dip 
dp 



VI' 



2pK' c 



(IV.12) 



The solution of this equation, going through the point 
(0, 0) on the phase plot, is 



-tp 



K 

V" 



(IV.13) 



In Fig. El we have plotted the phase curves obtained 
numerically for a toy model with the Lagrangian density 
p = \pp 2 — \ {p 2 + 2). As one can see from Fig. 
the parabolalike phase curve p s , given by the formula 
l|IV.13Jl . is the separatrix going through the fixed-point 
solution p(t) = 0. Moreover, this figure confirms that 
there are no phase curves intersecting the critical line by 
finite (p. 
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differentiable function, then for \(p\ <C 1 we have 



Figure 5: The typical behavior of the phase curves in the 
neighborhood of the critical line where K (ip) = (here p axis 
) is plotted for the case when K' c > and VI < 0. Horizontal 
dashed lines are the analytically obtained separatrices p>± and 
(0, u±) are the points of transition. 




Figure 6: The typical behavior of the phase curves in the 
neighborhood of the critical line where K (ip) = (here p axis 
) is plotted for the case when K c > 0, V c = 0, and V c " > 0. 



B. General differentiable K(ip) 

The models we are going to discuss below belong to 
the more general class of models for which the function 
K(ip) has zero of an odd order 2n + 1 (where n > 0) 
at p = 0. For the linear in X Lagrangians, we have 
e,xx = 0; therefore, if n > then K' c = and e° x = 0. 



That is why the general analysis made in the Sec. IIII CI 
does not work for this case. Therefore it is interesting to 
investigate on this simple example whether the desired 
transition could be possible for models not covered by 
our former analysis. If K(tp) is a sufficient many times 



K(<p) 



K 



(2n+l) 



(2n+ 1)! 



P' 



2n+l 



(IV.14) 



where Kc is the (2n + l)th ip derivative of K at 
ip = 0. If there is a phase curve, crossing the critical line 
at a finite nonvanishing ip c , then integrating both sides 
of the equation of motion III V. 41 we obtain 



Pc - Pi 



H(ip)dtp - 



(IV.15) 



Here (<fin,<pin) is a point on the phase curve in the neigh- 
borhood of the critical line. The first integral on the 
right-hand side of Eq. (|IV.15I) is always finite, whereas, 
as follows from the relation ljIV.14Jl . the second integral 
is definitely divergent, if e c v ^ 0. This divergence con- 
tradicts to our initial assumption: tp c - finite. There- 
fore we again obtain the condition 



l|IIL22jl . which re- 
stricts the possible intersection points on the critical line 
in the sense that in the other points, where the condi- 
tion does not hold, not only the second derivative (p does 
not exist, but there are no solutions pit) at all. More- 
over, it is clear that the condition ljIV.7Jl is not enough 
for the existence of the solutions intersecting the critical 
line. Thus, if the order of V'(ip) exceeds the order on ip 
of K'((p) for \(p\ <g; 1, then one can neglect V'(<p) and 
the integral IjlV.lFiJl has the logarithmic divergence (note 
that we do not consider the points ip c = because as 
we already know the transition does not occur via these 
points). When the order of V'(p) is lower than K'(<p) 
(and therefore lower than the order of K(ip)), we can ne- 
glect ip 2 K'(p) and the integral IHV.15JI has a power-low 
divergence. Finally, if the functions K'{p) and V'((p) 
have the same order on p for \p\ <g; 1 and are of opposite 
signs in a sufficient small neighborhood of p> = 0, then 
one can find an appropriate finite value ip^^Q for which 
the divergence on the right-hand side of Eq. ljIV.15Jl is 
canceled. One would expect that at this point the phase 
curves intersect the "critical" line and the dark energy 
changes the sign of (w + 1). Below, we give the direct 
calculation of these <p c and the phase curves in a neigh- 
borhood of them. Suppose that the order of the functions 
(V(p) — V c ) and K(ip) is (2n + 1) and there exist their 
derivatives of the order (2n + 2). Then for the p deriva- 
tive of the energy density we have in the neighborhood 
of the supposed intersection point (0, p c )- 



1 p 2n 

2 (2n)! 

P 



(2n- 



(p 2 c K c 2n+1) + K (2 " +1) ) + 2tpcK c 2n+1) 8Lp 

i) 0^< 



(2n+2) 



(2n+2; 



(IV.16) 



whereas the denominator pK(p) of the second integral on 
the right-hand side of Eq. iHV.lfiJl has the order (2n + 1) 
on p. The only possibility to get rid of the divergence 
in the integral under consideration is to assume that the 
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first term in the brackets in the asymptotic ijlV.lfit for 
e. v is zero. Therefore the possible crossing points are 
given by 



<Pc 



u± = ±\ 



(2n+l) 



(2n+l) ' 



K, 



(IV.17) 



Taking into account only the leading order on ip and Sip 
in the denominator and the numerator of Eq. I|IV.5|I , we 
obtain differential equation for the phase curves in the 
neighborhoods of the intersection points (0,u±); 



d ^=A ± -(2n+l) 6 A (IV.18) 

dp ip 



where 

A± = -3H C 



j^{2n+2)y{2n+l) _ ^(2n+l) y(2n+2) 



K, 



(2n+l) 



u± 



(IV.19) 

The solutions of this equation are given by the formula 



Sip 



2n + 2' 



2 »+ x = const. 



(IV.20) 



which is a generalization of formula IIIV.9II . Similarly to 
the case n — (K' c > 0) the solutions, on which the tran- 
sition occurs, have the measure zero in the phase curves 
set. Therefore we infer that the dynamical transition 
from the phantom states with w < — 1 to the usual with 
w > — 1 (or vice versa) is impossible. 

Now we would like to mention the models, for which 
V'(tp) is one order higher on <p than K'Up) for small ip. 
From the asymptotic expression for £ v ijlV.lfijl and the 
relation, giving the possible values of <p c ljIV.1511 . we see 
that the only point on the critical line which could be 
reached in a finite time is ip c = 0. Therefore, as we have 
seen in Sec. lIII Al the transition is impossible. The phase 
curve going trough the fixed-point solution ip(t) = is a 
parabola given by the generalization of Eq. l|IV.13ll : 



-ip 



K, 



(2n+l) 



(2n+2) 



Vr 



(IV.21) 



If V'((p) is more than one order higher on ip than K'((p), 
then as we have already mentioned ip c = and the tran- 
sition is impossible as well. 



C. Pole-like K(<p) 

In this subsection, we briefly consider the case when 
the function K(tp) has a pole of an odd order, so 
K — (/j~ 2n_1 , where n > 0, for \tp\ < 1. This kind of 
functions K(p) is often discussed in the literature in 
connection with the k— essence models (see [10]). Let 
us keep the same notation as in subsection IIV Al The 



potential V{p) can not have a pole at the point ip = 0, 
because, if it were the case, either the energy density e 
or the pressure p would be infinite on the critical line. 
In order to obtain finite values of the energy density 
e and pressure p, it is necessary to assume that the 
system intersects the critical line at ip = 0. But, as 
we have already seen in Sec. IIII Al the dark energy 
cannot change the the sign of (w + 1) at the points = 0. 

Thus, we have shown that in the particular case of 
the theories described by the linear in X Lagrangians 
p(ip,X) = K(tp)X — V(cp), which are differentiable in 
the neighborhood of (K(ip) and V(cp) differentiable 
but not necessary linearizable) the results, obtained for 
linearizable functions £,x,£,(p, and p t x , hold as well. This 
gives rise to hope that the same statement is true for the 
general nonlinear in X Lagrangians as well. Especially 
we have proven that, if the construction of the linear in 
X Lagrangian allows the transition, then the transitions 
always realize on a pair of the phase curves. One phase 
curve corresponds to the transition from w > — 1 to w < 
— 1 while another one realizes the inverse transition. This 
pair of phase curves obviously has the measure zero in the 
set of trajectories of the system. Therefore we infer that 
the considered transition is physically implausible in this 
case. 



V. SCALAR DARK ENERGY IN OPEN AND 
CLOSED UNIVERSES IN THE PRESENCE OF 
OTHER FORMS OF MATTER 

In the previous sections, we have seen that the desir- 
able transition from w > — 1 to w < - lis either impos- 
sible or dynamically unstable in the case when the scalar 
dark energy is a dominating source of gravity in the flat 
Friedmann universe. Let us now investigate whether this 
statement is true in the presence of other forms of matter 
and in the cases when the Friedmann universe has open 
and closed topology. 

Following Ref. [19], the effective sound speed c s is given 
by the same Eq. l|IL14ll for the flat, open, and closed uni- 
verses. Therefore, if the dark energy is the dominating 
source of gravitation (in particular this means that the 
energy density of the dark energy e ^ ) , then the anal- 
ysis made in Sec. IIII Bl is applicable to open and closed 
universes as well as to the flat universe. 

If the dark energy under consideration interacts 
with ordinary forms of matter only through indirect 
gravitational-strength couplings, then the equation of 
motion JII.1.3II can be written in the following form: 



<P£,X + "ipHp^x + E,tp = 0, 



(V.l) 



where merely the Hubble parameter depends on the spa- 
tial curvature and other forms of matter. This depen- 
dence is given by the Friedmann equation: 



H 2 + 



1 



(V.2) 
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where Si is the total energy density. It is obvious that 
the points on the plot (<p, ip) considered in the most of 
this paper do not define the whole dynamics of the sys- 
tem anymore and therefore do not define the states of 



the whole system. The analysis made in Sees. 
IIII CI and IIVI leans only on the behavior of the scalar 
filed <p and its first t derivative ip in the neighborhoods 
of their selected values, namely, such as where some of 
the conditions p.x — 0, ip = 0, or e t x = etc. hold. 
For these conditions, the contributions into the equation 
of motion ijV.lJl coming from the other forms of mat- 
ter and spatial curvature would be of a higher order and 
therefore are not important for the local behavior of ip 
and the problem as a whole. In fact, the value of the 
Hubble parameter did not change the qualitative futures 
of the phase curves considered in Sees. IIII AlIIlTH and 
IIVI To illustrate this statement, we plot the trajectories 

\2 



of the system p = ^<ptp 2 — \ ((<p — l) -4 + |J (it is the 

same system that we considered in previous subsection) 
in presence of dust matter for various values of the initial 
energy densities of the dust (see Fig. 01 . The only thing 
that is important is that H ^ 0. The universe should not 
change the expansion to the collapse and the plot of the 
scale factor a(t) should not have a cusp directly at the 
time of the transition. Thus, we infer that the most of 
our analysis is applicable to a more general physical situ- 
ation of a Friedmann universe filled with various kinds of 
usual matter, which interact with the dark energy only 
through indirect gravitational-strength couplings. More- 
over, if the interaction between the dark energy field ip 
and other fields does not include coupling to the deriva- 
tives V M y, then the obtained result holds as well. 



VI. CONCLUSIONS AND DISCUSSION 

In this paper we have found that the transitions from 
W > — 1 to w < — l(or vice versa) of the dark energy 
described by a general scalar-field Lagrangian p(tp, V^y) 
are either unstable with respect to the cosmological per- 
turbations or realized on the trajectories of the measure 
zero. If the dark energy dominates in the universe, this 
result is still robust in the presence of other energy com- 
ponents interacting with the dark energy through nonk- 
inetic couplings. In particular, we have shown that, un- 
der this assumption about interaction, the dark energy 
described by Lagrangians linear in (V M y>) cannot yield 
such transitions even if it is a subdominant source of 
gravitation. 

Let us now discuss the consequences of these results. If 
further observations confirm the evolution of the dark en- 
ergy dominating in the universe, from w > — 1 in the close 
past to w < — 1 to date, then it is impossible to explain 
this phenomenon by the classical dynamics given by an 
effective scalar-field Lagrangian p((p,V ^ip). In fact, the 
models which allow such transitions have been already 
proposed (see e.g., (HI H2> H3] anc ^ other models from the 





<j> fl 9 =o.mr 




Figure 7: Numerically obtained trajectories of the dark energy 
described by a Lagrangian linear in X are plotted for the cases 



and Q, v = O.lfV. 



Ref. |2J|) but they incorporate more complicated physics 
then the classical dynamics of a one scalar field. 

If observations reveal that w < — 1 now and if we dis- 
regard the possibility of the transitions, then the energy 
density of the dark energy should grow rapidly during 
the expansion of the universe and therefore the coinci- 
dence problem becomes even more difficult. Thus, from 
this point of view the transitions considered in this pa- 
per would be rather desirable for the history of the uni- 
verse. As we have shown, to explain the transition under 
the minimal assumptions of the nonkinetic interaction of 
dark energy and other matter one should suppose that 
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the dark energy was subdominating and described by a 
nonlinear in (V^) 2 Lagrangian. Thus, some nonlinear 
(or probably quantum) physics must be invoked to ex- 
plain the value w < — 1 in models with one scalar field. 

The second application of our analysis is the problem of 
the cosmological singularity. To obtain a bounce instead 
of collapse, the scalar field ip must change its equation of 
state to the phantom one before the bounce and should 
dominate in the universe at the moment of transition. 
Otherwise, if the scalar field was subdominant then it is 
still subdominant after the transition as well, because its 
energy density decreases during the collapse, while the 
other nonphantom forms of matter increase their energy 
densities. The disappearing energy density of <p does not 
affect the gravitational dynamics and therefore does not 



lead to the bounce. On the other hand, as we have al- 
ready proved, a dominant scalar field tp described by the 
action without kinetic couplings and higher derivatives 
cannot smoothly evolve to the phantom with w < — 1. 
Therefore we infer that a smooth bounce of the nonclosed 
Friedmann universe cannot be realized in this framework. 
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